Structural Studies of Decaying Fluid Turbulence: Effect of Initial Conditions 
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We present results from a systematic numerical study of structural properties of an unforced, 
incompressible, homogeneous, and isotropic three-dimensional turbulent fluid with an initial energy 
spectrum that develops a cascade of kinetic energy to large wavenumbers. The results are compared 
with those from a recently studied set of power-law initial energy spectra [C. Kalelkar and R. Pandit, 
Phys. Rev. E, 69, 046304 (2004)] which do not exhibit such a cascade. Differences are exhibited 
in plots of vorticity isosurfaces, the temporal evolution of the kinetic energy-dissipation rate, and 
the rates of production of the mean enstrophy along the principal axes of the strain-rate tensor. 
A crossover between non-'cascade-type' and 'cascade-type' behaviour is shown numerically for a 
specific set of initial energy spectra. 
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PACS numbers: 47.27.Gs 



I. INTRODUCTION 



Studies of structural properties (in space) of a tur- 
bulent fluid may lead to an understanding of dynami- 
cally significant processes in turbulence. Homogeneous 
and isotropic turbulence, within the context of statis- 
tically steady flows, has been the subject of extensive 
laboratory^ and numerical^, IS studies. These 
studies show that the turbulent fluid has a kinetic en- 
ergy spectrum with a wavenumber range that exhibits a 
power-law with exponent equal to —5/3, in accordance 
with a phenomenological theory due to Kolmogorov 
In the statistical steady state, visualization of the vortic- 
ity field reveals that regions of intense vorticity are orga- 
nized in slender filaments 0,0, IS El- The vorticity vector 
is observed to preferentially align |jj with the eigenvector 
associated with the intermediate eigenvalue of the strain- 
rate tensor of the fluid. By contrast, systematic spatial 
studies of decaying turbulence are scarce and restricted 
to numerical studies [SE1 of vorticity isosurfaces with ini- 
tial energy spectra of the type that develop a cascade of 
kinetic energy to large wavenumbers. 

A recent study 0, investigated the decay of unforced, 
incompressible, homogeneous, and isotropic three- 
dimensional magnetohydrodynamic turbulence from 
power-law initial conditions. The study was a general- 
ization of results [nj obtained for the corresponding fluid 
case. In particular, it was shown both analytically and 
numerically, that for the power-law initial energy spec- 
trum E(k,to) ~ k (k = |k| is the magnitude of the wave 
vector and t = in is the choice of virtual origin of time), 
the kinetic energy E(t) was found to decay as t' 1 and 
the integral length scale Lit) was found to grow as t 5 . 
Such power-law initial conditions are of interest in the 
astrophysical context of the decay of 'primordial' energy 
spectra 12J. It was shown numerically lOj, that such an 
initial energy spectrum does not develop a cascade to 



large wavenumbers; however, structural properties had 
not been probed. 

In this paper, we present results from a pseudospectral 
direct-numerical simulation (DNS) of the unforced, in- 
compressible, three-dimensional Navier-Stokes equations 
and contrast spatial results obtained from an initial en- 
ergy spectrum that develops a cascade (Section II), with 
results from the power-law initial spectrum (Section III) 
mentioned above. Plots of vorticity isosurfaces are found 
to differ, with distinct filaments in regions of intense vor- 
ticity, lacking in the case with power-law initial spectrum. 
The temporal evolution of the kinetic energy-dissipation 
rate and the rates of production of the mean enstrophy 
along the principal axes of the strain-rate tensor, are also 
found to differ. However, preferential alignment of the 
vorticity vector with the intermediate eigenvector of the 
strain-rate tensor is found to be unchanged. We also 
compare some classical results on mean enstrophy pro- 
duction in decaying turbulence. We show numerically 
that a crossover between the above-mentioned cases may 
be realised with an appropriate choice of initial condi- 
tions. 

The Navier-Stokes equations in vorticity form are 



~Dt 



= S, 



dxjdxj 



(1) 



where D/Dt = dt + v jdj is the material derivative, v 
is the kinematic viscosity, $y = l/2(djVi + diVj) is the 
strain-rate tensor, and u>i = tijkdjVk is the vorticity (e^-fe 
is the Levi-Civita tensor), i,j, k — 1,2,3, with a sum- 
mation implicit over repeated indices. We enforce the 
incompressibility condition diVi = 0. 



II. CASCADE- TYPE SPECTRUM 



A. Numerical Method 



'Electronic address: kalelkarQphysics. iisc.ernet. in 



We use a pseudospectral method^] to solve Eqs. (JIJ 
numerically, in a cubical box of side 2-7r with periodic 
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boundary conditions and 192 3 Fourier modes. We do 
not address issues pertaining to the scaling of higher- 
order structure functions of velocity differences or inves- 
tigate dissipation-scale properties, and believe that our 
spectral resolution is adequate for the types of studies 
that we have carried out. For the temporal evolution, we 
use an Adams-Bashforth scheme (step size St = 10~ 3 ) 
with double-precision arithmetic and set v = 1CP 5 . We 
include a hyperviscous term of the form i//jV 4 v in Eqs. 
JTJ, with v h — 10 -6 . We show explicitly in Section 11(C), 
that our results are unaffected by inclusion of hypervis- 
cosity. The initial velocity field is taken to be v(k, to) ~ 
k 2 e~ k e t6k , with 9^ random variables distributed uni- 
formly between and 2tt. This corresponds to an ini- 
tial energy spectrum E(k,t Q ) ~ fc 4 e~ 2fe " (with E(k,t) = 
|v(k, £)| 2 , the one-dimensional spectrum), which is a con- 
venient choice that develops a cascade to large wavenum- 
bers (see below). We measure time in units of the ini- 
tial large-eddy turnover time tq = 2ir/v® ms (here To 
equals 4.23), v° rms = [(£k |v(k, t )| 2 )] 1/2 is the root- 
mean-square value of the initial velocity, with the di- 
mensionless time r = £/tq (t is the product of the num- 
ber of steps and 8t). We define Reo = 2Trv^. ms /v to be 
the value of the initial 'box-size' Reynolds number (here 
Reo equals 932965). Our results are obtained for times 
to < t << t„, where is the time at which the (growing) 
integral scale L(t) = ((£ k |v(k, t)\ 2 /k)/ £ k |v(k, t)f) 
becomes of the order of the linear size of the simulation 
box. For times t > t*, finite-size effects which might well 
be non-universal, modify the numerical results, and are 
not considered here. 

In Fig. ^ we show some preliminary results that serve 
as a check of our numerical method and parameter val- 
ues (which were chosen to ensure linear stability of the 
numerical scheme). Figure ^a) shows on a log-log plot, 
the scaled kinetic energy spectrum fc 5//3 £'(fc, r) as a func- 
tion of the wavenumber k. On starting with the spec- 
trum specified above, a cascade of energy is seen to large 
wavenumbers. The plots are equispaced in time with a 
temporal separation of r = 0.24. The plot with open cir- 
cles is calculated at cascade completion at dimensionless 
time t = t c = 0.71, and shows a wavenumber range (for 
1 ^ k ^ 10) that exhibits the well-known —5/3 power- 
law 8] . Upon cascade completion, the shape of the en- 
ergy spectrum does not change appreciably (except at 
large wavenumbers where it falls), but the kinetic en- 
ergy decays. In Fig. ^b), we plot the normalized ki- 
netic energy E(t)/E as a function of the dimensionless 
time r, which is found to decay monotonically[^. In Fig. 
\H,c), we plot the normalized kinetic energy-dissipation 
rate e(r)/eo [e(t) = J2kk 2 \v(k,t)\ 2 } as a function of the 
dimensionless time t. The kinetic energy-dissipation rate 
peaks0,0| at r = r c , corresponding to cascade comple- 
tion in the energy spectrum, and decreases thereafter. 
The turbulence may be considered as 'fully developed' at 
t = r c and our spatial results (see below) will be calcu- 
lated at this instant of time. 




FIG. 1: (a) Log-log plot of the temporal evolution of the 
scaled kinetic energy spectrum k 5 ^ 3 E(k,r) as a function of 
the wavenumber k at temporal separations of r = 0.24, with 
E(k, to) ~ k 4 e~ 2k . The plot with open circles is calculated 
at cascade completion, at dimensionless time r = t c = 0.71. 

(b) Plot of the normalized kinetic energy E{t)/Eq as a func- 
tion of the dimensionless time r, with E(k,to) ~ k i e~ 2k . 

(c) Plot of the normalized kinetic energy-dissipation rate 
e(r)/eo as a function of the dimensionless time r, with 
E(k,t ) ~ fc 4 e- 2fe2 . 



B. Vorticity 

The dynamics of vortex structures in a turbulent flow- 
field is governed by the response of the vorticity vector to 
the strain-rate tensor. In Fig. [3(a) , we plot the normal- 
ized probability distribution P(|w|) of |a->| at cascade com- 
pletion. The distribution is found to peak at \uj\ — 1.5er 
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(cr denotes the standard deviation) and has an expo- 
nential tail for \ui\ > 1.5a with a fit P(\u\) ~ e"^ 1 ", 
P = 0.05 ± O.Of , a = 1 .00 ± 0.01 (error-bars from a 
least-squares fit). In Fig. Efb), we plot iso-|w| surfaces 
for the isovalue \u>\ = (\u>\) (angular brackets denote a 
volume average) at cascade completion, which appear 
to be crumpled sheet-like structures (found throughout 
the isovalue range [{\u)\) — a, {\u>\) + cr]). At early times 
r << r c , regions of intense vorticity (with the isovalue 
| a; | = (|a>|) + 2a) were found to be sheet-like (see Fig. 
121 a)). It is plausible that a small volume fraction^] of 
such structures, subject to shear and stretching, 'roll-up' 
to form the filaments visualised in Fig. EJb) at cascade 
completion. The filaments are found to be of diame- 
ter of the order of the grid spacing, and have a contour 
length that occasionally extends nearly to the linear size 
of the simulation box[15| . We choose to quote dimensions 
of the structures relative to the (fixed) box-size and the 
grid-spacing, since both the dissipative and the integral 
length scales vary in time, in decaying turbulence. Iso-|o;| 
surfaces in the range |w| < ((H) — a) were not found to 
exhibit any particular structure at cascade completion. 
The above results are in accordance with earlier stud- 
ies of statistically steady 0, 0, 0, H and decaying^, |j| 
turbulence. 



C. Strain- Rate Tensor 

At each grid point, we compute the eigenvalues Ai, 
A2 and A3 (with the convention Ai > A2 > A3) of the 
strain-rate tensor Sy, and the corresponding orthonor- 
mal eigenvectors ei, e2, and e^. In a constant-density 
flow, incompressibility requires that Aj = 0, with 
Ai > 0, A3 < 0, the sign of the intermediate eigen- 
value A2 being indeterminate. In Fig. Ufa), we pl°t the 
normalized probability distribution -P(A,) of the eigen- 
values Ai at cascade completion. We find that A2 has 
a positive mean. The statistically preferred ratio of the 
mean strain-rates (Ai) : (A2) : (A3) was found to equal 
4.9 : 1 : —5.9 at cascade completion [T^. 

In Fig. Efb), we plot the normalized probability dis- 
tribution P[cos(uj, of cosine of the angle between u> 
and the eigenvectors e», at cascade completion. We ob- 
serve that uj is preferentially aligned (or antialigned) 
with eigenvector ei (a peak-to- valley ratio ~ 5.1 was 
found) corresponding to eigenvalue A2 of minimum rela- 
tive mean strain-rate (disregarding the sign in the mean 
ratio, which merely indicates an extensional or compres- 
sional strain on the fluid element). The vorticity is also 
found to be preferentially perpendicular to eigenvector e% 
corresponding to the principal compressive eigenvalue A3 , 
while no specific angular relationship is observed relative 
to eigenvector e\. We could confirm an earlier result Q 
that the preferential alignment of w with ei is already 
pronounced during the cascade process, much prior to 
the appearance of distinct filaments in regions of intense 
vorticity. In Fig. Ufc), we plot P[cos(w, ej)] at cascade 
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FIG. 2: (a) Semilog plot of the normalized probability distri- 
bution P(M) of |w| (w is the vorticity) at cascade completion, 
with E(k, to) ~ k 4 e~ 2k . The dashed-line plot is a normalized 
Gaussian distribution for comparison. 

(b) Plot of iso-|(j| surfaces for the isovalue \lu\ — {\uj\} at cas- 
cade completion, with E(k,to) ~ fc 4 e~ 2fe . 

completion, for a run with the same parameter values 
as in Section 11(A), but without hyperviscosity (here 
v = 10~ 4 ). We find that the inclusion of hyperviscosity 
introduces no qualitative changes (a peak-to-valley ratio 
of 5.3 was found), and we believe that effects on struc- 
tural properties are minimal. The positive mean of A2, 
the ratio of the mean strain-rates, and the preferential 
alignment of u> with e2 is consistent with correspond- 
ing observations in statistically steadyQ and decaying^ 
turbulence. 



D. Mean Enstrophy Production Rate 

In a classical work[T^. Betchov derived an equation 
for the rate of production of the mean enstrophy (to 2 ) 
in unforced, incompressible, homogeneous, and isotropic 
turbulence, viz. 

^~^-i(GI&)> < 2 > 

where the angular brackets denote a volume average. 
From Eq. it is evident that production of (cu 2 ) re- 
quires a predominantly negative value of A1A2A3 in the 



4 




FIG. 3: (a) Plot of iso-|cj| surfaces for the isovaluc \lu\ = 
{\lu\} + 2<t (<t is the standard deviation) at dimensionless time 
t « T c , with E(k,t ) ~ fc 4 e" 2fe2 . 

(b) Plot of iso-|o;| surfaces for the isovalue \u\ = (\uj\) + 2a at 
cascade completion, with E(k,to) ~ fc 4 e~ 2ft . 



fluid. In Fig. [^a), we plot the normalized probability 
distribution P(AiA2As) of A1A2A3 at cascade completion 
and find that the distribution has a negative mean (with 
a skewness equal to —7.83), a result that has already 
been noted in an early numerical study [lflj. If A1A2A3 
is positive, we have one positive and two negative eigen- 
values: a vortex element is stretched along one direction 
and squashed along the other two, forming a slender fil- 
amentary structure. If A1A2A3 is negative, we have one 
negative and two positive eigenvalues: a vortex element 
is stretched along two directions and squashed along the 
third, forming a sheet-like structure. In homo gen eous 
turbulence (A1A2A3) is found to be negative 0, |2(|, in- 
dicating a predominance of sheet-like structures, which 
accords with Fig. |2Ib). In Fig. |3{b), we plot — (A1A2A3) 
as a function of the dimensionless time r. We find that 
the magnitude of — (A1A2A3) peaks at cascade comple- 
tion (r = 0.71), which is consistent with the dimension- 
less time at which the kinetic energy-dissipation rate is 
a maximum (cf. Fig. ^c)). In homogeneous turbulence, 
the term — (A1A2A3) can be shown [17j to equal the ex- 
pression {ujiUj-jSij) = ' e i) 2 )- In Fig. GJc), we 
plot (Xi(uj ■ ei) 2 }, the rates of production of the mean en- 
strophy along eigenvectors ei, and find that the produc- 
tion rate is largest along eigenvector ei corresponding to 





FIG. 4: (a) Plot of the normalized probability distribution 
P(Ai) of the eigenvalues \i of the strain-rate tensor Sij, 

i,j — 1,2,3, at cascade completion, with E(k,to) ~ k 4 e~ 2k . 

(b) Plot of the normalized probability distribution of cosine of 
the angle between u and the eigenvectors of Sij at cascade 
completion, with E(k, to) ~ k 4 e~ 2k . 

(c) Plot of the normalized probability distribution of cosine of 
the angle between u and the eigenvectors d at cascade com- 
pletion, with E(k,to) ~ k 4 e~ 2k and without hyperviscosity. 



eigenvalue A2 of minimum relative mean strain-rate (see 
above), a remarkable result that appears to be new. 
From Fig. I^b) we observe that the mean enstrophy pro- 
duction rate — (A1A2A3) has an upper bound. In the same 
work^^i Betchov invoked formal mathematical inequal- 
ities and incompressibility to derive the upper bound 

|(A 1 A 2 A3)|<^=((A 2 + A2 + A2) 3 / 2 ). (3) 
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FIG. 5: (a) Semilog plot of the normalized probability dis- 
tribution P(AiA2As) of A1A2A3 at cascade completion, with 

E(k,to) ~ fc 4 e _2fc . The dashed-line plot is a normalized 
Gaussian distribution for comparison. 

(b) Plot of — (A1A2A3} (see Eq. as a function of the di- 

mensionless time r, with E(k,to) ~ fc 4 e _2fc . 

(c) Plot of the rates of production of the mean enstrophy 
(Aj(td • ei) 2 ) along the eigenvectors et of S»j, as a function of 
the dimensionless time r, with E(k,to) ~ fc 4 e _2fc . 



The upper bound is the largest rate of production of 
mean enstrophy compatible with the requirements of 
isotropy, homogeneity, and incompressibility. In Fig. [ft] 
we plot inequality © 21j as a function of the dimension- 
less time r and find that the inequality is satisfied to 
a greater degree during the cascade process, and for di- 
mensionless times r > t c asymptotes to the value « 0.57. 
Dynamical effects suggest an inequality that is stronger 
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FIG. 6: Plot of Betchov's inequality B = 
3v^|(AiA 2 A3)|/((A 2 + A! + Al) 3/2 } < 1 as a function of 

the dimensionless time r, with E(k,to) ~ k 4 e~ 2k (open 
circles) and E(k,to) ~ k (asterisks). 



than the one due to Betchov. 



III. POWER-LAW SPECTRUM 
A. Numerical Method 

We use the numerical scheme, spectral resolution, step 
size, viscosities, and boundary conditions as specified in 
Section 11(A). The initial velocity field is taken to be 
v(k, to) ~ fc 1 / 2 e l ^ k , with (f>k random variables distributed 
uniformly between and 2w, which corresponds to the 
initial energy spectrum E(k,t ) ~ k (with E(k,t) = 
|v(k, t)\ 2 ). Here, the dimensionless time To equals 13.9 
and the initial 'box-size' Reynolds number Reo equals 
350181. Our results are obtained for times t < t « £*, 
and we choose to calculate our spatial results at dimen- 
sionless time t = t c = 0.7l[23 which is equal to the di- 
mensionless time at which the kinetic energy-dissipation 
rate peaks, on starting with a 'cascade- type' spectrum as 
shown in Fig. ^c). 

In Fig. EJa), we show on a log-log plot, the kinetic en- 
ergy spectrum E(k, r) as a function of the wavenumber 
k. The plots are equispaced in time with a temporal 
separation of r = 0.24 and the plot with open circles is 
calculated at r = t c . The spectrum does not cascade to 
large wavenumbers or exhibit a wavenumber range with 
a —5/3 power-law at any stage of temporal evolution. In 
Figs. EK D ) and (c), we plot the normalized kinetic energy 
E(t)/Eq and the normalized kinetic energy-dissipation 
rate e(r)/eo, as a function of the dimensionless time r 
(plotted on a linear scale in order to compare with Figs, 
^b) and (c)). The kinetic energy E(t) is found to decay 
as a power-law (on a log-log plot) with an exponent equal 
to —0.91 ± 0.04, with error-bars from a least-squares fit. 
The exponent is theoretically predicted |TT| to equal — 1, 
and we believe the discrepancy is due to the low spectral 
resolution of our DNS. The normalized kinetic energy- 
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FIG. 7: (a) Log-log plot of the temporal evolution of the ki- 
netic energy spectrum E(k, r) as a function of the wavenum- 
ber k at temporal seperations of r = 0.24, with E(k,to) ~ k. 
The plot with open circles is calculated at dimensionless time 
r = r c = 0.71. 

(b) Plot of the temporal evolution of the normalized kinetic 
energy E(t)/Eq as a function of the dimensionless time r, 
with E(k,to) ~ fc. 

(c) Plot of the normalized kinetic energy-dissipation rate 
e(r)/eo as a function of the dimensionless time r, with 
E(k,t ) ~ fc. 



FIG. 8: (a) Semilog plot of the normalized probability distri- 
bution P(|w|) of \u>\ (to is the vorticity) at dimensionless time 
t = t c , with E(k,to) ~ fc. The dashed-line plot is a normal- 
ized Gaussian distribution for comparison, 
(b) Plot of iso-|w| surfaces for the isovaluc \lu\ = at 
r = r c , with E(k,to) ~ fc. 



B. Vorticity 

In Fig. 02a), we plot the normalized probability dis- 
tribution P(|w|) of \oj\ (uj is the vorticity) at r = r c . 
The distribution peaks at |w| = 1.5<r as in Fig. HJa), 
however it does not exhibit a stretched-exponential tail. 
In Fig. Efb), we plot iso-|w| surfaces for the isovalue 
\oj\ = at r = r c , which appear to be shredded sheet- 
like structures (observed throughout the isovalue range 
[(\uj\) — a, (\uj\) + a]). At early times r << r c , regions of 
intense vorticity (with the isovalue |w| = (\lo\) + 2a) were 
found to be structure-less 'blobs' of length of the order 
of the grid spacing (see Fig. E{a)), in contrast to the 
sheet-like structures in Fig. I^a). At r = r c , isosurfaces 
of intense vorticity (in Fig. E^b)), appear to be roughly 
ellipsoidal in shape with semiaxes of the order of the grid 
spacing and differ markedly from the coherent filaments 
in the 'cascade-type' case shown in Fig. Efb). 



C. Strain- Rate Tensor 



dissipation rate e(r)/eo does not exhibit a peak (cf. Fig. 
Die)) and decays monotonically. 



In Fig. llUf a). we plot the normalized probability dis- 
tribution P(Xi) of the eigenvalues A; of the strain-rate 
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FIG. 9: (a) Plot of iso-|cj| surfaces for the isovaluc \u>\ = 
+2<j (cr is the standard deviation), at dimensionless time 
t « t c , with E(k,to) ~ ft. 

(b) Plot of iso-|w| surfaces for the isovalue = + 2<r at 
r = r c , with E(k, to) ~ fc. 



FIG. 10: (a) Plot of the normalized probability distribution 
P (Ai) of the eigenvalues Ai of the strain-rate tensor Sij, i,j = 
1, 2, 3, at dimensionless time r = r c , with E(k, to) ~ fc. 
(b) Plot of the normalized probability distribution of cosine of 
the angle between to and the eigenvectors e% of at r = t c , 
with E(k,t ) ~ k. 



tensor Sij at r = r c . We find that A 2 has a positive 
mean. The statistically preferred ratio of the mean eigen- 
values (Ai) : (A 2 ) : (A3) was found to equal 4.6 : 1 : —5.6, 
which differs only marginally from that obtained in the 
'cascade-type' case. In Fig. HOf b). we plot the normalized 
probability distribution of cosine of the angle between u> 
and the eigenvectors e, of Sij at r = r c . As in Fig. 1Kb), 
we observe that u> is preferentially parallel (or antiparal- 
lel) to 62 (with a reduced peak-to- valley ratio ss 4.1) and 
perpendicular to e3, while no specific angular relation- 
ship is observed relative to e\ (though we note a small 
probability of alignment at cos(w, ei)=±l). 



of the dimensionless time r, and find that the production 
rates decrease monotonically along all the principal axis 
directions, in remarkable contrast to Fig. OJc). In Fig. 
El we plot inequality JSJ as a function of the dimension- 
less time r, and find that the inequality monotonically 
approaches the same asymptotic value « 0.57 as in the 
'cascade-type' case. 



E. Crossover 



D. Mean Enstrophy Production Rate 

In Fig. Illf a). we plot the normalized probability dis- 
tribution P(\\\i\z) of A1A2A3 at r = t c - As in Fig. 
|S{ a), we find that the distribution has a negative mean 
(a relatively smaller skewness equal to —3.90 was found). 
In Fig. ITlT b). we plot the value of — (A1A2A3) (see Eq. 
(J5J) as a function of the dimensionless time r and find 
that it decreases monotonically, in contrast to Fig. EJb). 
In Fig. Hir e), we plot (Aj(w • e^) 2 ) , the rates of production 
of the mean enstrophy along eigenvectors as a function 



From Figs, ^c) and[7{c), it is evident that some rep- 
resentative set of initial energy spectra may be utilised to 
exhibit a crossover from non-'cascade-type' to 'cascade- 
type' behaviour, with the kinetic energy-dissipation rate 
serving as a diagnostic. In particular, one may choose 
the set E(k,to) ~ k q e~ k with parameter q. In Fig. 1121 
we plot e(r)/eo as a function of the dimensionless time 
r (with To calculated using the initial energy spectrum 
with q = 0.4). On varying q in the range [0.4,0.8], we 
find that e(r)/eo exhibits a smooth crossover at q ~ 0.6 
from non-'cascade-type' to 'cascade-type' behaviour. 
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FIG. 11: (a) Semilog plot of the normalized probability distri- 
bution P(AiA2Aa) of A1A2A3, with E(k,to) ~ k. The dashed- 
line plot is a normalized Gaussian distribution for comparison. 

(b) Plot of — (A1A2A3} (see Eq. as a function of the di- 
mensionless time r, with E(k,to) ~ k. 

(c) Plot of the rates of production of the mean enstrophy 
(Ai(td • ei) 2 } along the eigenvectors a of Sij, as a function of 
the dimensionless time r, with E(k,to) ~ fc. 
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FIG. 12: Plot of the normalized kinetic energy-dissipation 
rate e(r)/eo as a function of the dimensionless time r, with 
E(k,t ) ~ k q e- k ", 0.4 < q < 0.8. 

isosurfaces, the temporal evolution of the kinetic energy- 
dissipation rate, and the rates of production of the mean 
enstrophy along the principal axes of the strain-rate ten- 
sor. However, preferential alignment of the vorticity vec- 
tor with the intermediate eigenvector of the strain-rate 
tensor is found to be unchanged. A crossover between 
non-'cascade-type' and 'cascade-type' behaviour is shown 
numerically for a specific set of initial energy spectra. 
Our study shows that decaying turbulence is strikingly 
different for different types of initial conditions. Initial 
conditions that lead to a cascade yield features similar to 
those obtained in the statistically steady case. However, 
power-law initial conditions that do not lead to an energy 
cascade exhibit several qualitatively new features which 
we have described above. 
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IV. CONCLUSIONS 



To summarize, we have presented results from a sys- 
tematic numerical study of structural properties of an un- 
forced, incompressible, homogeneous, and isotropic tur- 
bulent fluid with an initial energy spectrum that devel- 
ops a cascade to large wavenumbers. The results are 
contrasted with those from the power-law initial spec- 
tra investigated in Refs. 0, [nj, which do not exhibit 
such a cascade. Differences are noted in plots of vorticity 
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